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The single-site dynamical mean field theory approximation to the double exchange model is found 
to exhibit a previously unnoticed instability, in which a well-defined ground state which is stable 
against small perturbations is found to be unstable to large-amplitude but purely local fluctuations. 
The instability is shown to arise either from phase separation or, in a narrow parameter regime, 
from the presence of a competing phase. The instability is therefore suggested as a computationally 
inexpensive means of locating regimes of parameter space in which phase separation occurs. 

PACS numbers: 71.10.+W, 71.27. +a, 75.10.-b, 78.20-e 
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I. INTRODUCTION 



Dynamical mean field theory (DMFT) has been widely 
applied to many strongly correlated electron systems pj. 
Since DMFT takes local quantum fluctuations into ac- 
count, it is especially successful for models whose many 
body effect comes from the on-site interaction, like 
HubbardQ or Kondo^Q (double exchange) model. 
Correlated systems often exhibit different phases which 
are quite close in energy, and this proximity can lead to 
phase separation, which is often important for electronic 
physics [5|- Phase separation is in principle a "global" 
property of the phase diagram and requiring substan- 
tial effort to establish: one must compute the free en- 
ergy over a wide parameter range, and then perform a 
Maxwell construction. In this paper we show that within 
the single-site DMFT formalism a straightforward calcu- 
lation at a fixed parameter value can reveal the presence 
of phase separation. Specifically, we find that at zero 
temperature, the DMFT can give a ground state which 
is stable against small perturbations but is unstable to a 
large amplitude local perturbation; at non-zero temper- 
ature the standard methods simply fail to converge to a 
stable solution. By computing the free energy and per- 
forming a Maxwell construction we show that for wide 
parameter ranges this instability occurs in the regions 
in which phase separation exists. In a narrow parame- 
ter regime it signals instead the onset of a different, but 
apparently uniform, phase. We therefore propose that 
the instability of the DMFT equations can be used as an 
approximate, computationally convenient estimator for 
the boundaries of the regimes in which phase separation 
occurs. 

The balance of this paper is organized as follows. We 
first present the model and then a zero temperature dy- 
namical mean field analysis explicitly showing the insta- 
bility. In section III we calculate the full T=0 phase dia- 
gram in the energy-density plane, establish the regime of 
phase separation via the usual Maxwell construction, and 
extend the treatment to T > 0. In section IV we discuss 
the implications of this instability. Finally in section V 
we present a brief conclusion. 



II. DOUBLE EXCHANGE MODEL AND 
DYNAMICAL MEAN FIELD APPROXIMATION 

A. Double Exchange Model 

In this paper, we consider the single orbital "double 
exchange" or Kondo lattice model of carriers hopping be- 
tween sites on a lattice and coupled to an array of spins. 
This model has been studied by many authors and con- 
tains important aspects of the physics of the " colossal" 
magnetoresistance manganites and is also solvable in a 
variety of approximations, permitting detailed examina- 
tion of its behavior. Here we use it to investigate the 
physical meaning of a previously unnoticed instability of 
the dynamical mean field equations. 

The model is defined by the Hamiltonian 

H = 4 «( c £<t c .j> + h - c -) +^JS l - c+ a a a p Ci t0 

i 7 Ot.fi 



(1) 



with i,j labeling the sites and the Si denoting the spins. 



We assume the spins are classical ( 



Si , 5*7 



0) and are 



of fixed length. We choose the convention = 1. 

The hopping ty defines an energy dispersion Sk and 
thus a density of states D{e) = J d d k/ {2-K) d 8(e — In 
our actual computations we specialize to the d — » oo limit 
of the Bethe lattice, for which D{e) = V^t 2 - e 2 /(27rf 2 ) 
because the availability of convenient analytical expres- 
sions allows us to accurately compute the small difference 
between free energies of different states. We choose en- 
ergy units such that t = 1. 

We also note that the ground state properties of the 
model may be straightforwardly obtained, because for 
any fixed configuration of the spins the model is quadratic 
in the fermions and easily diagonalizable. 



B. Dymanical Mean Field Method 

We now present the dynamical mean field analysis 
of this model. In the single site dynamical mean field 
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method 1], one neglects the momentum dependence of 
the self energy. The properties of the model may then 
be calculated by solving an auxiliary quantum impurity 
model, along with a self consistency condition. The quan- 
tum impurity model corresponding to EqnQ is specified 
by the partition function 



Zimp — / dS e" 



(2) 



with A = Trlog [a^ai ~ J cos9(a^ — ajj — J 2 ] where the 

trace is over frequency, and cos 9 = z ■ S. S is deter- 
mined by a spin-dependent mean field function a(u). In 
a magnetic phase, a-f ^ aj.. Note that the assumption 
of classical core spins means that J dS denotes a sim- 
ple scalar integral over directions of the core spin £>, and 
that no Berry phase term occurs in the argument of the 
exponential. 

The Green function Gimp and self energy S(lj) of the 
impurity model are given by 



5 log Z. 



imp 



5a a (u)) 



= a a - G 



imp,u 



(3) 



a is fixed by requiring the impurity Green's function Gi mp 
equals to the local Green's function of the lattice prob- 
lem. The form of the self consistency equation depends 
on the state which is studied. For a ferromagnetic (FM) 
state, it is 

G im p,a = + /i - £j - Scr^))" 1 (4) 

fccBZ 

while for a 2 sublattice antiferromagnetic (AF) state, 

-l 



G 



fccRBZ 



(5) 



where £1 (2) = iu> n + /i — Sf m and the k -sum is over tje 
reduced Brillouin zone (RBZ). The two equations become 
equivalent in the paramagnetic (PM) state where Ef = 
S| 

The solution of Eqs [3] to [5] determines the magnetic 
phase, the single particle properties, and the free energy. 
In particular, in the dynamical mean field approximation 
the Gibb's free energy is 00 



N 



Qimp — TTrlagG a (iw n ) - TTi log G^ 1 (iuj n , k) (6) 



where the trace is over the spin and lattice degree of 
freedom. The Helmholtz free energy is F(n) = + fj, N. 
At zero temperature, the ground state energy (Helmholtz 
free energy) is 



E 
N 



Tr[G(iu n ,k)(iLLi n + fj,)] 



(7) 



We also note that the solution a a (iuj n ) defines an ef- 
fective potential for the core spin, which depends on the 
angle 9 between the core spin and local magnetization 
direction, so that 



J imp 



D -E eff (cos9)/T 



with 

E eff (cost 



-T^Tlog 



J cos 9(a^ — oj.) + J 2 



(8) 



(9) 



C. Phase Boundaries and Maxwell Construction 

The model is known to exhibit ferromagnetic, spiral 8] 
and commensurate antiferromagnetic phases. For our 
subsequent analysis, an accurate determination of phase 
boundaries will be important. We therefore present here 
a few calculational details. 

We require the T = phase boundary separating the 
ferromagnetic and spiral phases 4J 8] . The energy of a 
spiral state may most easily be found by performing a 
site-dependent spin rotation to a basis in which the spin 
quantization axis is parallel to the local spin orientation. 
The problem may then be easily diagonalized. For the 
infinite dimensional Bethe lattice one finds, for a diagonal 
spiral of pitch </>, that in the rotated basis, the local Green 
function is given byQ 

G~V) = w + fi - a J - t 2 cos 2 (^)G a - t 2 sin 2 (^)G s 

(10) 

where is the angle between two nearest neighbor mag- 
netization, and tilde is used for the spiral states. 

To locate the T = FM/spiral second order phase 
boundary, it suffices to expand ground state energy 
E(fi,(f)) (EqnQ) to second order in (p. The energy dif- 
ference between FM and spiral states is 



dE\ 



on 



N 



—d(j)=[-\ T^2iLo n dG(iuj n ,fi) (11) 



with dG equaling to 



d G a — G 



G a 



t 2 (Go — Go) 
t 2 - Go 2 



(12) 



where a = —a. The FM/spiral phase boundary is deter- 
mined by dE\N — 0. 

The model also exhibits phase separation in some 
regimes. To determine the boundaries of the regime 
where phase separation occurs, we use the DMFT 
method to compute the Helmholtz free energy as a func- 
tion of occupation number F(n) (Eqn©) and then per- 
form the Maxwell construction. An example is shown in 
Fig ^ We find that in fact that over much of the phase 
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FIG. 1: Energy as function of particle density at J = 1, 
T = 0, illustrating Maxwell construction. Heavy solid line: 
FM state. Light solid line: 2 sublattice AF state. Dashed 
line: Maxwell construction interpolatiion. 
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(a) FM Instability 
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(b) AF Instability 



FIG. 2: (a) Effective potential for the core spin (Eqn©) for a 
ferromagnetic state calculated at T — 0, J = 8, and densities 
shown. As we increase the occupation number from 0.88 to 
0.90, the minimum of E e / /(cos 6) changes from cos 9 = — 1 to 
1 which indicates the FM DMFT solution becomes unstable 
when n > 0.89. (b) Effective potential for the core spin for 
an anti-ferromagnetic state calculated at T = 0, J = 1, and 
densities shown. At n = 1, the minimum of E e // happens 
at cos# = —1, while at n = 0.995 the minimum changes to 
cos9 = 1. The AF DMFT solution only exists at n = 1. 



diagram a phase separation between FM and n = 1 AF 
states preempts the formation of spiral or n < 1 AF 
state. The general structure of our phase diagram agrees 
with earlier work|g| [To| . but the precise locations of phase 
boundaries differ by roughly 10%. 



where < ... > means the angular average with respect 
to the weight function e -E?s}( c ° sS )/ T ) w ith E eff defined 
in Eqn© . At zero temperature, the only contribution 
of the angular average is from the absolute minimum 
of E ef f (cos 0). To find the DMFT solution at T = 0, 
one first assumes the absolute minimum of E e ff (cos 0) 
occurs at a fixed value, for example cos(9) = —1, ob- 
tains GT = Oj,(f) — (+)J from Ean (|13|l . and gets 
the self energy £f(n = +(— )J from Eqn®. Finally, 
one uses the a obtained by the above procedure to cal- 
culate E e /f (cos 9) to see if the minimum is located at 
the point originally assumed. Note that different ground 
states (FM, AF..) enter the above procedure only via the 
self consistent equation Eqn|@} (or EqnJSJ). 

Fig El shows that as density is increased at fixed large 
J, the self consistency breaks down, in an unusual man- 
ner: E e ff (cos 9) remains locally stable (slope around the 
assumed minimum cos(9) = —1 remains positive) but 
the global minimum of E e /f moves to cos(9) = 1. In 
the regime where this phenomenon occurs, no solution of 
the DMFT equations exists. Any initial solution we have 
considered leads to a similar inconsistency (as is shown 
in panel b of Fi§J2for the case of anitferromagnetism) . 

This instability is also manifest at T > 0. As T is 
decreased at fixed n, the convergence becomes slower and 
below some temperature T*(n), no stable solution can 
be found for a J-dependent range of n. The absence of a 
solution for some range of n can be seen in a different way 
by solving the model as a function of chemical potential fi 
Fig|S{a) shows that as \x is increasedat fixed low T, a first 
order transition occurs to a paramagnetic state, with a 
corresponding jump in n. Associated with the first order 
transition is a coexisting region in which two solutions are 
locally stable (FM with lower n and PM with higher n); 
the DMFT equations correspondingly have two solutions, 
which one is found depends on the initial seed. The solid 
and dashed lines in Fig |3fa) are obtained from initial 
seeds close to FM and PM states respectively. 

The absence of convergence may be understood from 
the density dependent effective potential, shown e.g. in 
Fig E{b). One sees that as n is increased, E(cos9 = 1) 
decreases; this is a precursor of the effect shown in Fig 
Efa)- Indeed, the curve E(cos9) is reduced by an n- 
dependent scale factor. For n larger than a critical value 
(here n ~ 0.905), E(cos9 = 1) is small enough relative 
to the temperature that this region begins to contribute 
to < cos 9 > , lowering the maximum m that can be sus- 
tained and destabilizing the ferromagnetic solution. 



III. DMFT INSTABILITY 



IV. INTERPRETATION 



In this section we show that the DMFT equations ex- 
hibit an apparently previously unnoticed instability. We 
begin with T = 0. From Eqn©, the Gimp is 



G 



imp,t(l) 



=< 



H(T) ~ (+)^cos( 



a-fdj — Jcos6'(a| — ajj — J 2 



> 



(13) 



We argue in this section that the DMFT instability 
documented in the previous section is a manifestation 
of competing instabilities (primarily phase separation) 
in the original model. To establish this we show in the 
panel (a) of FigEla T = phase diagram in the density- 
coupling plane. The dash-dot line shows the phase sep- 
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FIG. 3: (a) n(p) at J = 8 and T = 0.01, showing two phase 
behavior. Solid line: obtained from the initial seed close to 
the FM state. Dashed line: obtained from the initial seed 
close to the PM state, (b) E ef f(cos9) at J = 8, T = 0.01, 
n = 0.88(solid line), 0.89 (dot line), and 0.90(dashed line). 
The circle and dagger represent the scaled effective potential 
for n = 0.88 and 0.89 which are almost indistinguishable to 
the E e ff at n = 0.90. 



FIG. 5: (a) Phase Boundaries for J=8 in n — T plane. Solid 
line: the phase separation boundary obtained by Maxwell 
construction. Dotted line: Curie temperature. Dashed line: 
the boundary of DMFT instability, (b) Maxwell construction 
for J — 8, T — 0.008 - the difference between free energy and 
interpolating line. The ground state is FM when n < 0.86, 
and FM/AF phase separaion for 0.86 < n < 1. 




FIG. 4: (a) T = phase boundaries. Solid line: DMFT in- 
stability above which the DMFT solution becomes unstable 
Dot line: spiral instability line. The system is spiral(FM) 
above(below) this line. Dashed line: AF/FM phase separa- 
tion, (b) Expansion of the region 0.8t < J < 2t. For n, J 
in the triangle bounded below by the solid line and above by 
the dashed and dot lines, the ground state is not known. 



aration boundary obtained from the global energy com- 
putation; for n above this line the model phase separates 
into an n = 1 AF and an n < 1 FM state. The dotted 
line shows the phase boundary between uniform FM and 
spiral states. Finally, the heavy solid line shows the re- 
gion above which the FM DMFT solution is unstable at 
T = 0. When J is large enough that the FM state is 
fully polarized (J > 2t), we see that the DMFT insta- 
bility line follows the phase separation line, but is inside 
the region of phase separation. We therefore suggest that 
in this region the DMFT instability is a consequence of 
phase separation and this DMFT instability line can be 
used as a rough estimate of the real phase separation 
boundary. 

When J < 0.8*, the DMFT instability indicates the 



presence of a spiral state with lower energy than the fer- 
romagnetic state. For 0.8* < J < 1.7* (FigHfb)), there 
exists a narrow region of n where none of the uniform 
phases we considered solve the DMFT equations and the 
Maxwell constructions seem not to indecate phase sep- 
aration. We believe that in this region there exists a 
uniform non FM/AF/spiral/paramagnetic state (either 
the ground state or the phase separation beteen FM and 
that state) which we do not know yet. 

At T > the situation is similar. The DMFT instabil- 
ity is contained inside the regime of phase separation. For 
example, we show in Fig^a) the phase diagram and the 
range of DMFT instability in the density-temperature 
plane for J = 8*. The heavy line shows the boundary 
of the regime of phase separation obtained by Maxwell 
construction: for < T < 0.011, the phase separation 
is between FM and AF(n = 1); for 0.012 < T < 0.02, 
the phase separation is between PM and AF(n = 1); 
for 0.011 < T < 0.012, the phase separation is either 
PM-AF(n = 1) or FM-PM[ll| according to the location 
n, T relative to the homogenous Curie temperature (dot- 
ted line). The dashed line shows the region where the 
DMFT solution fails to converge at that given density 
n (the DMFT equation has stable solution for all fi, see 
Fig EDO)- For n > 0.95, the DMFT instability line de- 
notes the temperature below which (a) the paramagnetic 
state is linearly unstable to antiferromagnetic and (b) no 
stable antiferromagnetic solution exists (except n = 1). 

FigEJb) shows the results of a Maxwell construction 
for J = 8 and that at T = 0.008, presented as the differ- 
ence between calculated free energy F(n) and the inter- 
polating line I(n) = F{n = 1) + F{n "\~_ F n ^ =1) {I - n) 
with n* — 0.86. Phase separation is seen to occur 
for n* < n < 1, while the DMFT instability range is 
0.895 < n < 0.95 and 0.96 < n < 1. 
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V. CONCLUSION 

We have found an instability in the ferromagnetic 
DMFT equation for the single site double exchange 
model and shown that this instability corresponds to the 
FM/AF phase separation when the coupling J is larger 
than half bandwidth (2t) and to another ground state 
(spiral) in the small coupling region. There exists a 
small window, around intermediate J, where no stable 
FM DMFT solutions exist while the spiral or phase sep- 
aration is not the ground state, and we believe there is 
a non FM/AF/Spiral/Para ground state existing in this 



region. We have presented evidence that the instability is 
a signal, obtained from a calculation at a fixed parameter 
value, of the existence of an instability (typically phase 
separation) which normally is established via a global 
computation, comparing free energies at many different 
parameter values. We therefore propose that the DMFT 
instability is a computationally convenient way to esti- 
mate the boundary of phase separation. 

We thank Dr. Satoshi Okamoto for many helpful dis- 
cussions. This work is supported by DOE ER46169 and 
Columbia University MRSEC. 
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